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Abstract 

Populations are seldom completely isolated from their environment. Individuals in a particular 
i-Sh ' geographic or social region may be considered a distinct network due to strong local ties, but 

will also interact with individuals in other networks. We study the susceptible-infected-recovered 
(SIR) process on interconnected network systems, and find two distinct regimes. In strongly- 
coupled network systems, epidemics occur simultaneously across the entire system at a critical 
infection strength /3 C , below which the disease does not spread. In contrast, in weakly-coupled 
network systems, a mixed phase exists below /3 C of the coupled network system, where an epidemic 
occurs in one network but does not spread to the coupled network. We derive an expression for the 
network and disease parameters that allow this mixed phase and verify it numerically. Public health 
implications of communities comprising these two classes of network systems are also mentioned. 

PACS numbers: 64.60.aq, 87.10.Mn, 89.75.Da 
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I. INTRODUCTION 



Complex network models of the interactions in human society have been used to under- 
stand many problems in epidemiology [l-sj. These models have generally assumed that all 
of the nodes interact on a single network with a single degree distribution. Even when these 
degree distributions allow for large heterogeneities — as in the case of scale- free networks 
where hubs with large numbers of connections can arise — the assumption remains that every 
node is part of a single network and is represented by a single underlying topology. In re- 
ality, however, societies are composed of many interconnected networks, as in Fig. [H which 
may be communities within a larger population or separate systems entirely. A disease can 
spread through the network of direct personal contacts, via the water utilities network, and 
through travel from city to city over highway or airline networks. These interconnected net- 
work systems may be comprised of different types of nodes, which may have degrees drawn 
from distinct degree distributions, and may have different connectivities between them. Real 
world examples of these systems can be seen in a 2009 study by Stehle et. al. which found 
a three fold difference in interaction time between students inside and outside of their own 
class 10 1 . Other studies have shown similar patterns [ill . 11^ . Human- animal interact- 



ing network systems are also of great importance. The H5N1 variant of influenza spreads 
through the network of birds, and from them to individuals in the network of humans that 
work or live closely with them. While no current mechanism exists for efficient spreading 
from human to human, there is substantial concern that such a mechanism will evolve 13 ) . 
Human-mosquito-human transmission for the P. knowlesi malaria strain is also a worrying 
concern [l^ . 

Interconnected network systems have been of interest to researchers in numerous different 
ways flo1-[l8|. Interconnected dependency networks, where failure in nodes in one network 
causes failures of dependent nodes in the other network, exhibit failure cascades, where the 
cross network dependencies result in a network much more easily fragmented than single net- 
works of the same degree distribution jisl . Interconnected power networks, where transport 
capacity and failure vulnerability are competing properties, were examined and an optimal 



level of interconnection found 



20|. Networks without dependencies, such as interconnected 



social networks, where populations exist at city, state, and national levels, have also been 
examined. In these networks, the level of movement between cities (the interconnections be- 
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FIG. 1: An interconnected network system with two networks: A and B. Nodes have intranetwork 
links within their own network, but also internetwork links connecting them to the other network. 



ween them) have been shown to affect the epidemic transition on the metapopulation level 
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22j, although in this case the low-level networks were treated in a mean field fashion, 



classified only by rate equations and infection numbers, with no internal network features. 
In addition, the percolation threshold in interacting networks was found to be lower than in 



23]. 



single networks, with a giant cluster appearing for a smaller total number of links 

In this work we consider two interconnected networks (or, alternately, interconnected 
communities within a single, larger network). We pose the question: Under what conditions 
will an epidemic spread only on the sub-networks, with minimal isolated infections on other 
network components, and under what conditions will it spread across the entire intercon- 
nected network system? Depending on the parameters of the individual networks and their 
interconnections, connecting one network to another can have a profound or a small effect 
on the spread of an epidemic. Identifying the conditions in which these cases occur is vital 
to our understanding and management of epidemic processes. 

We define two different interconnected network regimes, strongly and weakly coupled, 
and find the interaction strength value separating these two regimes. Our primary result 
is to show that in the strongly-coupled case, we find that all networks are simultaneously 
either disease free or part of an epidemic, while in the weakly-coupled case a new "mixed" 
phase can exist. In this mixed phase, the disease is epidemic on only one network, and not in 
other networks, despite the interconnections. The applications to public health are straight- 
forward. If two neighboring communities comprise a strongly-coupled network system, then 
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an outbreak in any community is cause for immediate concern in the other. Due to this, 
in the strongly-coupled case it becomes important to pursue a strategy of communication 
and joint action between public health agencies, and perhaps even intervention from a single 
agency with higher authority 



II. MODEL 



In this section, we consider the case of only two interconnected networks of equal size, 
but it is easily possible to extend the model to an arbitrary number of networks of any size. 
We form our interconnected network systems in the following way: 

1. Generate two networks, A and B, with their own intranetwork (A -H- A and B -B- B) 
degree distributions Pa(^) an d -fs(^) according to the standard Molloy-Reed config- 



uration model 



24]. 



2. Draw a degree from the internetwork (A -B- B) degree distribution, Pab^), for all the 
nodes in both networks. 

3. If the total degree assigned to nodes in network A is not equal to the total degree 
assigned to nodes in network B, randomly reassign a node in B until the total numbers 
in each network are equal. 

4. Randomly connect nodes in network A to nodes in network B to form the intercon- 
nected network system. 

This method generates random, uncorrelated, interconnected network systems with specified 
inter- and intra-network degree distributions. While this method works for any arbitrary 
degree distribution, P x (k), we present results only for random Poissonian degree distribu- 
tions. 

The susceptible-infected-recovered (SIR) epidemic model is used here to study the effects 
of interconnected network structure on epidemic threshold. The SIR model is well estab- 
lished and describes diseases such as HPV, seasonal influenza, or H1N1 2l|, |25|]. In this 
model, each node has three possible states: susceptible (s), infected (i), or recovered (r). 
Each node begins in state s, except for a single node in one network chosen to be in state 
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i. Nodes in state i infect their neighbors in state s with probability (3 at each time step, 
changing them to i. Nodes enter state r after spending a recovery time t r in state i. 

In order to find the threshold for an epidemic, we can think of epidemic spreading as a 
bond percolation process [26, 127]] on a network. In bond percolation, links between nodes 
are activated with a certain probability p. If this probability is greater than a certain critical 
value, p c , then a giant cluster emerges, where the existence of a path between any two nodes 
is almost certain. In a disease-spreading model, nodes infect their neighbors, "activating" 
the links between them with a certain probability, and a disease reaches nodes through this 
entire network above a certain critical value, /3 C , just as in the case for percolation. 

In complex networks, this critical threshold for percolation if all potential links are acti- 
vated is k = 2. Here k is the expected number of nearest neighbors that a node chosen by 
following an arbitrary link will have, and is calculated from the ratio between the second and 
the first moments of the degree distribution: k = (k 2 ) / (k). For k > 2, a giant cluster exists, 
while for k < 2 only small isolated clusters exist. If some subset of bonds is activated at 
random with probability p, a giant cluster appears at a critical value of p c = — 1) (28]. 

The SIR model likewise has an epidemic phase transition at a critical (3 = (3 C below which 
the disease remains confined to the local neighborhood of the initial infection, and above 
which the disease spreads throughout the network. This transition from the disease-free 
phase to the epidemic phase depends on the average number of secondary infections per 
infected node becoming larger than one. This allows the long-term survival of the disease, 
as the infection density will grow over time on average, and thus ensure that the epidemic 
spreads to a large fraction of the population. In our problem, the expected number of 
susceptible neighbors that a node has when it just becomes infected is given by n — 1, since 
the total expected number of neighbors is k, and one of them must be excluded as the infected 
parent from which the current node descended. The transmissiblity = 1 — (1 — j3) tr is the 
probability to infect a neighbor before recovery. The mean number of secondary infections 
per infected node is thus Nj — (ac — l)Tp. The infection will die out if each infected node 
does not infect on average at least one replacement so, for a very large network, the critical 
point is given by the relation (k — l)Tp = 1. The single network model exhibits only a single 
transition at (3 C given by Q| 

&(«) = !- [I-Ck-I)" 1 ] 17 *-. (1) 
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In the interconnected network model, the behavior is more complicated, as the disease can 
potentially cause an epidemic in different combinations of the networks. The disease can 
either be in the epidemic phase in both networks, in the disease-free phase in both networks, 
or active in one network while the other remains disease free, called here the mixed phase. 
The boundaries of these phases are controlled by k a , Kb, and kt, where k a and Kb are 
calculated over the individual A and B networks, disregarding internetwork connections, 
and kt is calculated over the entire coupled network system, including intra- and inter- 
network links. 



III. STRONGLY-COUPLED NETWORK SYSTEMS 

We consider an interconnected network system to be strongly coupled if kt is larger than 
ka-i and Kb- For random networks, k = (k) + 1, and thus we may write kt in terms of the 
average degrees (&b), and (/cab) as follows: 

«T = [(k A ) 2 + (k A ) + (k B ) 2 + (k B ) + 2(k AB ) 2 + 2(k AB ) 
+2(k A )(kAB)+2(k B )(kAB)} 

[(kA) + (k B ) + 2(k AB )]- 1 . (2) 

Without loss of generality, we define network B as the more intraconnected network ((ks) > 
(kA))- For fixed network parameters (/c^) and we can then derive the critical interaction 
strength (/cab)c that separates strongly-coupled (n T > n B ) from weakly-coupled (k t < k b ) 
networks: 

{ k AB)c = • (3) 

In strongly-coupled network systems, we expect any epidemic to emerge simultaneously on 
networks A and B. Using Eq. [TJfor each of the three k, it can be seen that for the strongly- 
coupled case /3 c (kt), the critical value of (3 for the disease to emerge on the giant component 
formed by the entire interconnected network, is smaller than both (3 c (k a ) and /3 c (kb), the 
critical values of /3 for epidemics to spread on networks A or B ignoring internetwork links. 
As such, any pathogen virulent enough to spread in network A or B alone will have already 
caused an epidemic occurring across the interconnected network system. For this case, 
the disease spreads across the interconnected network system as a single network, with 
the internetwork connections bringing an epidemic into existence before any intranetwork 
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connections can do so independently; the mixed phase will not be seen. To support this, we 
plot the ratio of the largest connected infected cluster formed solely from nodes connected 
with intranetwork links, compared to the size of the largest connected cluster formed by 
nodes connected with all links, in Fig. |5J For a strongly-coupled network system, the relative 
size of the largest connected infected component contained entirely in a single network 
decreases initially, showing that the epidemic is occurring across the interconnected network 
system, not locally in one of the networks. Thus in the strongly-coupled case, epidemic 
spreading is enhanced due to internetwork connections, with epidemics occurring for less 
virulent diseases than would spread on either network alone (lower C .) 

One note is that for networks of identical intranetwork degree ((&a) — (^b)) (^ab)c = 0. 
That is to say, identical networks always form strongly-coupled network systems. This is 
in agreement with our findings that the phase diagram of strongly-coupled network systems 
is similar to that of single networks. An interacting network system formed by attaching 
the labels 'A' and 'B' to different halves of a single network would be such an example 
system, and one should not expect that this relabeling could have any effect on the physical 
properties, such as phase transitions, of that network. 

IV. WEAKLY-COUPLED NETWORK SYSTEMS 

If two networks are connected with (/cab) below the threshold value from Eq. EJ i.e. 
Kb > ^Ti we define the interconnected network system to be weakly coupled. From Eq. [2] 
this also gives n T > Turning again to Eq. [T]we see that /3 c (n B ) < /3 c ( k t) < Pc^a)- 

Epidemic spreading is a non-competitive process. Adding more links to a network can 
only increase the spread of an epidemic, never decrease it, as the chance of a node infecting 
its neighbors is constant regardless of degree. Thus, a disease with (3 above the individual 
epidemic threshold of network B ((3 c (kb)) will enter the epidemic phase on that network, 
regardless of the other network and the values of kt and k,a- If P is below (3 c (kt), however, 
the disease cannot spread to more than isolated small clusters of network A. Thus, in the 
weakly-coupled case, we expect to see a mixed phase, with the boundaries dependent on the 
values of /3 and (A;ab)- A mixed phase indicates that the addition of the interconnections 
between the two networks is only affecting epidemic spreading on the network with weaker 
intranetwork connections, with the epidemic on the network with stronger intranetwork 
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FIG. 2: Epidemics in strongly-coupled network systems spread across all networks, remaining 
confined to one network in weakly-coupled network systems. We plot Sb, the size of the largest 
connected cluster solely in network B, divided by St, the size of the infected cluster across the 
interconnected network system, for both strongly- and weakly-coupled network systems. The B 
only cluster decreases in relative size until criticality (/3 c (kb) = -048), showing that the epidemic 
spreads throughout both networks rather than remaining confined in B in the strongly-coupled 
system. By contrast, in the weakly-coupled system, the relative size of the B only cluster grows 
until j3 = P c {kt) = .054, showing that growth is localized in the more strongly coupled network. 
For the strongly-coupled network, (kX) = 1-5, (ks) = 2.5, and (/cab) = 2.5. For the weakly-coupled 
network, (k^) = 1.5, (ks) = 4.55, and {Hab) = 0.3. In all cases, Na = Nb = 10 4 and t r = 5. 

connections unchanged by the internetwork links. The weakly-coupled case in Fig. [2] shows 
this, with the largest connected cluster contained entirely in B becoming larger compared 
to the size of the giant component with increasing (3 until {3 c {k>t) is reached, indicating that 
the disease does not spread through connected regions of network A. 

If f3 is increased to above /3 c (kt), network B becomes capable of spreading the disease 
to network A, which now enters the epidemic phase, even for /3 < /3 c (ka)- This matches 
the work done by Leicht and D'Souza 20] where a giant cluster forms consisting of nodes 
in both networks, even when the less intraconnected network is below its own percolation 
threshold. We plot the full phase diagram for both weakly- and strongly-coupled networks 
in Fig. |3l showing the disease-free phase, the mixed phase, the epidemic phase, and the 
transition between weakly- and strongly coupled networks. The existence of this mixed phase 
is important in the real- world context of interacting networks, as the communities or systems 
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FIG. 3: The mixed phase disappears in the transition from weakly- to strongly-coupled network 
systems. Sample phase diagram for interacting network systems with = 1.5 and (ks) = 6.0 as 
a function of infection strength /3 and internetwork degree {J^ab), showing the two critical j3 and 
(&ab)c- In the weakly coupled case, below /3 c (kb) no epidemic occurs. For /3 c (kb) < (3 < /3 c (kt), 
there exists a mixed phase, where an finite fraction of network B becomes infected, but network 
A has only small infected clusters. Above /3 c (kt), an epidemic occurs across the entire network in 
both the weakly- and strongly-coupled cases. For this diagram, N A = N B = 10 4 and t r = 5. 

that comprise the components are likely to be governed by different bodies. If two cities, for 
example, together form a weakly-coupled network system, the more highly connected city 
can more safely disregard the links to, and response of, the less highly connected city, as the 
spread of the epidemic will depend on local parameters only. 

We performed Monte-Carlo simulations to verify this result. First, Fig. HI shows infection 
densities at different (3, corresponding to a horizontal sweep across the phase diagram seen 
in Fig. [3] at (/cab) = 0.1. The epidemic spreading first occurs at i3 c (kb), where the disease 
enters the epidemic phase and spreads through network B, while the infection density in 
network A remains negligible. This mixed region, in agreement with our predictions, occurs 
in the region /3 c (kb) < (3 < /3 c («/r). In this regime, network A plays no role in the spreading 
of the infection on network B. Above (3 c (kt), we see that the infection density in network A 
begins to rise, showing that the entire interconnected network system is now in the epidemic 
phase, as predicted. 

For networks approaching the strongly-coupled regime from below, the mixed phase is 
expected to be small, and thus difficult to identify from graphs such as the one in Fig. HI 
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FIG. 4: In the mixed phase, the two networks have separate transition values. Ratio of total number 
of infected Nj in each network to the size of the giant cluster Nqc f° r two weakly coupled networks 
with (kA) = 1.5, (ks) = 6.0 and (&ab) — •!■ The respective epidemic thresholds calculated from 
Eq. (prj are /3 c (k_b) ~ 0.035, (3 c (kt) ~ 0.0425. The infection can be seen to become epidemic in 
network B well before it does in network A. The network of networks has Na = Nb = 10 4 and 

tr- 5 . 



We thus examine not only the infection densities, but also the survival probability P(t), 
which is the probability of an infection started from a single infected site being active at 
a time t. Equivalently and more accessible from public health records, the distribution of 
time spans of reported outbreaks can be used. At criticality, the probability of an infection 
started from a single infected site remaining active at a later time t is expected to scale 
as P(t) ~ t~ l 29[. FigJH shows the survival probabilities of the networks comprising an 
interconnected with = /3 c (/cb) for both the strongly- and weakly-coupled cases. In both 
cases, Network B exhibits the expected t~ l fall-off in survivability with time that is expected 
of a system at criticality, indicating that Network B is actually undergoing a phase transition 
at the disease-free/mixed phase line. In the weakly-coupled case, however, the survival 
probability in network A does not fall off as expected, due to infrequent and non-epidemic 
instances of infections from network B. The slope of the survival probability for network 
A thus cannot be used directly to confirm when it enters the epidemic phase. However, 
if both networks are participating in an epidemic, the disease should be active in both 
networks at each time step. We thus introduce the survival probability gap (inset of Fig. |3j), 
AP(t) = min [[Pg(t) — Pa{1)] /P{t)], or the minimum relative difference in the likelihoods 
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FIG. 5: Epidemics exhibit critical survival only in one network for weakly-coupled networks. 
Infection survival probabilities P(t) on the individual networks with internetwork connectivity 
(&ab) = .1 at /? = /3 c (kb)- The survival probability in network A (lower curve, with +), the less 
connected network, is a small fraction of that in network B (upper curve with 0)> the more con- 
nected network. The survival probability in network B falls off as f , as expected of a system at 
criticality. Network A does not show a smooth decrease towards typical of a network much below 
criticality. The inset shows the relative difference in the survival probabilities, [Pb(£) — P J 4(i)]/P(t). 
The arrow indicates the minimum difference between the two curves, after the initial increase. Net- 
work parameters are (/ca) = 1-5, (fee) = 6.0, Na = Nb = 10 4 , and t r = 5. 

that each network will have any infected members (nodes in class i) present at time t. We 
use this quantity to measure the deviations of the survival probability in network A from 
the value that is obtained from a network at criticality (network B) and thus of how far 
away network A is from it's own critical point. 

In Fig. |6] we plot this survival gap at different /3, equivalent to vertical slices across the 
phase diagram seen in Fig. |3j We see that when (/cab) or P is increased to move outside 
the expected mixed phase region and into the epidemic phase, AP(t) goes to zero. In other 
words, at the mixed/epidemic phase line, network A is behaving identically to a network 
known to be at criticality, and thus can be said to itself be critical along that line. This 
confirms the assertion that there can only be a gap in survival probability when one network 
is in the epidemic phase and the other is not, i.e. in the mixed phase. Thus a non-zero 
survival probability gap can serve as a good predictor for the presence of the mixed phase. 

Lastly, we addressed the question of universality under different values of inter- and intra- 
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FIG. 6: Survival probability gap shows mixed phase boundaries. Fractional size of the minimum 
survival probability gap, AP(t)/P(t r ), (minimum distance between the two curves in Fig. [5] after 
the time t r has passed) between two interacting ER networks with (/ca) = 1-5 and (&b) = 6.0 at 
various infection strengths. From top to bottom j3 = /3 c (kb) = .0358, .038, .04, .042,/3 c (kt) = .044. 
The gap shrinks with increasing interaction and with increasing (3, matching Fig. [3j At /3 c (kt), 
the survival probability is the same in both networks for all but (&ab) = -01, where it remains 
distinct due to finite size effects. For all systems, Na = Nb = 10 and t r = 5. 

network degree, finding that along the disease-free/mixed phase transition line, the behavior 
of networks with different k is universal under appropriate scaling. FigJT] shows three different 
networks, all with > (fc^) and /3 = (3 c (kb)- Rescaling the survival probabilities by P(t r ) 
and plotting vs. k,t/^b instead of kt directly, the curves collapse, showing AP(t), and thus 
the mixed phase, disappear uniformly as the network approaches the strongly-coupled regime 
(kt > Kb and (/cab) > (^ab)c-) Near the critical point, AP(t) ~ [(/vr — kb)/kb}~ 1 (fit not 
shown.) This identical behavior implies survival probabilities in networks could be used 
as as a measure of network connectivities near criticality, as the latter may be difficult to 
obtain for social and biological networks, whereas information on the duration of an epidemic 
outbreak in various communities (from which P(t) can be estimated) is likely to be recorded. 
In addition, the survival probability gap persists well beyond {Ivab) = 1, for appropriate (/c^) 
and (ks)- For the system with (/c^) = 3.0 and (ks) = 12.0, (Jzab) c ~ 2.47. Even when every 
node in network A is connected to two or more nodes in network B, there can still be an 
order of magnitude difference at minimum in the likelihood of finding the disease active in 
the two networks; the mixed phase region is not confined to small (/c^b) only. 
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FIG. 7: Survival probability gap scaling is universal when rescaled by kb for different networks. 
The data collapses onto a single curve, showing this universal behavior. The ratio of kt to Kg 
determines AP(i) along the disease-free / mixed phase border. For the network systems with 
(ks) = 6.0, /3 = .0358, for (ks) = 12.0, j3 = .01725. Both are /3 c (kb) of the respective systems. 
Again, for all systems Na = Nb = 10 4 and t r = 5. 

V. CONCLUSIONS 

In summary, we introduced two classes for interconnected network systems, strongly cou- 
pled and weakly coupled, and studied the behavior of epidemics on them. In strongly-coupled 
network systems, epidemics occur always across the entire interacting network system , with 
the presence of interconnections enhancing epidemic spreading. In weakly-coupled network 
systems, a mixed phase exists where epidemics do not always occur across the full inter- 
connected network system, and interconnections affect only epidemic spreading across the 
less intraconnected network. We demonstrated the boundaries and behavior of the mixed 
phase numerically as well as analytically. Proper analysis of which groups of communities 
comprise strongly- or weakly-coupled systems could inform public policy and highlight the 
necessity of cooperation between different governing bodies or provide information about 
the epidemic danger of increasing interaction between human and animal populations. 
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